Vorticity production and survival in viscous and magnetized cosmologies 
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We study the role of viscosity and the effects of a magnetic field on a rotating, self-gravitating fluid, 
using Newtonian theory and adopting the ideal magnetohydrodynamic approximation. Our results 
confirm that viscosity can generate vorticity in inhomogeneous environments, while the magnetic 
tension can produce vorticity even in the absence of fluid pressure and density gradients. Linearizing 
our equations around an Einstein-de Sitter cosmology, we find that viscosity adds to the diluting 
effect of the universal expansion. Typically, however, the dissipative viscous effects are confined to 
relatively small scales. We also identify the characteristic length bellow which the viscous dissipation 
is strong and beyond which viscosity is essentially negligible. In contrast, magnetism seems to favor 
cosmic rotation. The magnetic presence is found to slow down the standard decay-rate of linear 
vortices, thus leading to universes with more residual rotation than generally anticipated. 
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I. INTRODUCTION 

Most, if not all, astrophysical systems rotate and it is 
conceivable that vorticity can play an important role in 
a number of situations. So far, however, the majority of 
the related studies has been confined to relatively small- 
scale systems. Although the question of universal rota- 
tion and how large this may be has been posed more than 
once (see, e.g., [l|-|6[ and references therein), in cosmol- 
ogy the role of vorticity is usually bypassed. One of the 
reasons for doing so is that typical inflationary scenar- 
ios are irrotational,since the vacuum fluctuations do not 
generate vorticity p}- Another, is that linear rotational 
perturbations are known to decay as the inverse square 
of the cosmological scale factor, which implies strong di- 
lution for any vortex-like distortions that could happen 
to exist. Dissipative effects during the subsequent radia- 
tion era, like Silk damping and neutrino free-streaming, 
may also have a negative impact on large-scale rotation, 
though (to the best of our knowledge) studies of this na- 
ture are not yet available in the literature. 

In cosmology, where the terms vorticity and rotation 
are used interchangeably, vorticity is usually generated 
from nonlinear effects. This is one additional reason 
for not taking vorticity into account in many cosmologi- 
cal perturbation studies. Inhomogeneous non-barotropic 
(e.g., non-adiabatic) media, for example,can generate ro- 
tation at the seco nd p erturbative level [H [t| , and so can 
fluid viscosity |10l - [l3 |. Nevertheless, resorting to nonlin- 
ear effects is not always necessary, since there are agents 
that can induce rotation still at the linear level. For in- 
stance, a coherent magnetic field ('-B-field') can act as a 
source of cosmic vorticity [HI, [TJ] . The opposite is also 
theoretically possible. Thus, following [la ], a number of 
authors have explored the likelihood that rotational dis- 
tortions can lead to cosmic magnetogenesis (see, e.g., fl6l — 
[HI and references therein). In addition, both the fluid 



viscosity and the -B-field can affect the evolution of rota- 
tional perturbations and thus change the standard pic- 
ture, which usually assumes non-magnetized, perfect me- 
dia. Therefore, in principle at least, the effects of the 
aforementioned sources could lead to more residual vor- 
ticity than previously anticipated. 

The present paper uses Newtonian gravity to study the 
effects of viscosity and those of a magnetic field on ro- 
tation. We adopt the covariant, Lagrangian approach to 
fluid dynamics, which has been proved helpful in several 
(mainly relativistic) cosmological studies (see [2(| for an 
introduction to the Newtonian version of the formalism). 
First, we look into the effects of fluid viscosity and then 
consider those of a coherent magnetic field. Our results, 
confirm those of earlier studies claiming that both of the 
aforementioned agents can act as sources of cosmic ro- 
tation. We show, in particular, that viscosity can only 
generate vorticity in media with an inhomogeneous den- 
sity distribution. This is not the case for the B-field, 
which does not (in principle at least) require the pres- 
ence of density perturbations. Assuming the Newtonian 
analogue of the Einstein-de Sitter universe as our back- 
ground model, we find that viscosity can only generate 
vorticity at second order. Magnetism, on the other hand, 
can do the same at the linear perturbative level. More- 
over, the irreducible decomposition of variables and equa- 
tions that the covariant formalism introduces, allows us 
to identify the agents leading to vorticity generation, as 
well as the way they do so. Viscosity, for example, sources 
rotational distortions via the coupling of the shear with 
certain types of density perturbations. Alternatively, the 
magnetic field acts primarily through the tension compo- 
nent of the Lorentz force. 

Our next step is to investigate how viscous and/or 
magnetized media can affect the linear evolution of ro- 
tational perturbations. In the first case, we find that 
viscosity always leads to vorticity dissipation, though the 
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effects are generalfy scale dependent. We identify, in par- 
ticular, a characteristic length beyond which viscosity is 
essentially negligible and below which its role becomes 
progressively more important. In particular, well in- 
side this "viscous scale" vorticity dissipates very quickly, 
while far beyond that scale the fluid kinematics are dom- 
inated by the background Hubble expansion. There, one 
recovers the standard inverse-square decay law for ro- 
tational perturbations. The size of the aforementioned 
viscous scale depends on the specifics of the cosmic fluid, 
such as the average velocity and the mean free-path of 
the species involved. When dealing with non-relativistic 
matter, however, the typical domain affected by viscosity 
is very small (in cosmological terms). 

In contrast to viscosity, the presence of magnetic fields 
seems to have a favorable effect on linear rotational per- 
turbations. In addition to generating vorticity at the 
linear order, the B-field slows down the standard de- 
cay rate of rotation in perturbed Friedmann-Robertson- 
Walker (FRW) universes. As mentioned earlier, vorticity 
typically depletes as a~ 2 , with a representing the cosmo- 
logical scale factor. When the linear magnetic effects are 
accounted for, however, the vorticity decay-rate changes 
to a~ 3 / 2 . The implication of this effect, which comes en- 
tirely from the magnetic tension, is that magnetized uni- 
verses, can have more residual vorticity than their non- 
magnetic counterparts. An analogous, conclusion was 
also reached indirectly, through the relativistic analysis 
of linear density inhomogeneities in magnetized cosmolo- 
gies [1J]. To the best of our knowledge, this is the first 
direct indication that magnetism "favors" both the gen- 
eration and the survival of cosmic rotation. 



II. NEWTONIAN SELF-GRAVITATING MEDIA 

A. Newtonian kinematics 

Consider a self-gravitating Newtonian medium and a 
velocity field v a , tangent to the flow lines of the fluid. 
The local kinematics is analyzed by splitting the gradient 
dbv a into its irreducible components according to [20] 

dbV a = - QS ab + (T ab + £abcU C , (1) 

where O = d a v a is the volume scalar, <j ao = d^v^ rep- 
resents the shear tensor and uj a = — cvslv a /2 defines the 
vorticity vector. 1 Also, S a b is the Kronecker symbol and 
s a bc is the Levi-Civita tensor. By construction, the vol- 
ume scalar determines the expansion or contraction of 
a fluid element (when positive or negative, respectively) 
and defines a representative length scale (a) by means 



of a/a = 0/3. The shear, on the other hand, describes 
changes in the element's shape (under constant volume) 
and u! a defines the local rotational axis. 

The above defined variables obey a set of three evolu- 
tion and three constraint equations that fully determine 
the fluid kinematics. In particular, the propagation for- 
mulae are @, [H| 

Q = -\Q 2 - 1 -np + d a A a -2{a 2 -lo 2 ) , (2) 



&ab = — g Qo ab - E ab + d( b A a ) - (T c ( a <7 C b ) - UJ( a U> b ) (3) 

and 

2 1 

w» = - - 6w a - - curL4 a + a ab uj , (4) 

where overdots denote convective derivatives (e.g. = 
d t Q + v b db&)- Here, p is the density of the matter (with 
k = 8irG being the rescaled gravitational constant), while 
2er 2 = a a b(J ab and 2uj 2 = uj a bio ab define the magnitudes of 
the shear and the vorticity, respectively. We also note the 
vector A a — v a +d a Q, which is the sum of the inertial and 
gravitational acceleration ($ is the gravitational poten- 
tial) and vanishes in the absence of non-inertial or non- 
gravitational forces. 2 Finally, the tensor E a b = d( a db)Q 
describes the tidal part of the gravitational field. 

Expressions ©-(HI) are supplemented by an equal 
number of constraint equations. These relate the gra- 
dients of the key kinematic variables and are given by 

d a uj a = , curler afc = -d( b u) a ) (5) 

and 

2 

d b <T a b = g <9 Q 6 + curlw a , (6) 

with curkxafc = £ c d{ a d c cr d b)by definition [ISIllI- The sets 
and ©-(O monitor the kinematics between two 
neighboring flow lines. To close the system one also needs 
the continuity equation, which reflects mass conservation 
and in our case is 

P = -e P) (7) 

together with the momentum-conservation formula. The 
exact form of the latter depends on the type (i.e. on the 
equation of state) of the matter and on the presence of 
additional sources. 



Round brackets denote symmetrization, square antisymmetriza- 
tion and angled ones indicate the symmetric and trace-free part 
of second-rank tensors (e.g. a a ^ = d^ v a) = d^ v a) ~ (0/3)<5 a (>)- 



2 Putting the inertial and the gravitational acceleration together 
means that there is no longer an explicit need for Poisson's equa- 
tion. The latter, however, has been employed when deriving the 
set of tSJ) — 13]) and takes the form d 2 <& = np/2 in our notation, 
with d 2 = d a d a representing the Laplacian operator. 
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B. Newtonian cosmologies 

Newtonian cosmology is expected to provide a good 
approximation of the post-recombination universe and 
on scales well inside the Hubble radius, where one can 
safely ignore the role of spacetime curvature. In New- 
tonian gravity, the relativistic FRW universes are repre- 
sented by models where the fundamental observers move 
along parabolic, elliptical and hyperbolic trajectories. 
These flow lines correspond to Euclidean (flat), spherical 
(closed) and hyperbolic (open) spatial geometry respec- 
tively. Given the homogeneity and isotropy of the Fried- 
mann universes, the only nonzero variables are scalars 
that depend solely on time. Vectors, tensors and spatial 
gradients vanish identically. Thus, the Newtonian ana- 
logue of a dust-dominated, fiat FRW cosmology (i.e. of 
the Einstein-de Sitter universe) is monitored by the set 



H 2 = \n Pl H 



-H z 



1 



and 



-3Hp. 



(8) 



(9) 



These relations are respectively known as the Friedmann, 
the Raychaudhuri and the continuity equations. Note 
that H = a/ a is the Hubble parameter, with a rep- 
resenting the cosmological scale factor. This implies 
that H — 0/3 in all FRW cosmologies. It is then rel- 
atively straightforward to show that (|5b) and © are 
the Einstein-de Sitter limits of © and respectively. 
Friedmann's equation, on the other hand, is simply an 
integral of (|5b). 

The system of © and © can be solved analytically 
with respect to time, leading to 



a oc t 



2/3 



H 



and up 
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(10) 



which offer an alternative description of the Einstein-de 
Sitter universe. In the following sections, we will con- 
sider linear rotational perturbations around such a back- 
ground, and study their evolution in the presence of fluid 
viscosity or a magnetic field. 



III. VISCOUS ROTATION 

A. Viscosity induced vorticity 

The rotational evolution of a self-gravitating Newto- 
nian fluid is governed by Eq. (j4]), which determines the 
twist between two neighboring flow lines. Note the first 
and the last terms on the right-hand side of that relation. 
These describe effects triggered by the relative motion of 
the fluid particles, analogous to those attributed to the 
Coriolis force [ll[. The second term, on the other hand, 
reflects the presence of non-inertial and non-gravitational 
forces (e.g. those triggered by pressure gradients, the 



presence of electromagnetic fields, etc) and is determined 
by the associated Navier-Stokes equation. In the case of 
a (non-magnetized) viscous fluid, the latter reads [2(| 



pA a = -d a p - d b ir a b 



(11) 



where p is the isotropic pressure and ir ab is the anisotropic 
pressure of the matter (with ir ab = n ba and ir a a = 0). 
Taking the "curl" of the above and then substituting 
into expression Q provides an evolution formula that 
monitors the rotational behavior of a Newtonian viscous 
medium. 

Fluid viscosity can be expressed in terms of kinematic 
(shear) viscosity via the phenomenological relation 



-2vpa ab , 



(12) 



with v representing the kinematic viscosity coefficient 
and v > [ll|. The latter immediately implies that 
v = n/(2pa), where 2n 2 — n ab ir ab defines the magnitude 
of the anisotropic pressure tensor. In addition, when the 
viscosity coefficient is constant, the divergence of Eq. (fT2")l 
ensures that 



d b ir ab = -2vo ab d b p - - vpd a Q - 2vpcurluj a , (13) 



where we have used the so-called shear constraint (see 
Eq. ([5]) earlier). 

Inserting the auxiliary relation (| 13[) into the right-hand 
side of (fTT|) , taking the "curl" of the resulting expression, 
using the constraint ([5^) and then substituting the out- 
come into Eq. (Q}, the latter gives 

2 1 

w 8 = - 5 ©w a - — e abc d b (In p) d c p 

3 2p 



ve abc d b [a cd d d (\np)]+ud 2 uj a 



(14) 



According to the above, in the absence of viscosity, the 
only way of generating vorticity is through the baroclinic 
(the second) term on the right-hand side of (fT4")l . The 
latter vanishes when p — p(p), namely in the case of a 
barotropic medium. The viscous nature of the fluid, on 
the other hand, can act as a source of rotation by ex- 
ploiting the coupling between the density and the shear 
gradients. Analogous conclusions have also been reached 
in [Tlj |. I n fact, after removing the baroclinic term and 
taking into account the difference in the definition of the 
vorticity vector, one can show that expression (|14p is 
identical to Eq. (1) in (Tlj . 

Let us consider the viscous source-term in the right- 
hand side of HU), t ne significance of which means that 
it deserves further attention. Recalling that curl<j b = 
£ cd(ad c<jd b) an d using constraint (JSJd), a straightforward 
calculation leads to the expression (see also Appendix [SJ 



e a bcd b [a cd d d (\np)] = 



- 2A b d {b uJ a) - A b e bcd d c a d a 
+ e abc X bd a c d , (15) 
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which decomposes the source term into its irreducible 
components. According to the above, the source term is 
not entirely rotation-free, implying that only a part of 
it can actually generate vorticity. Also, A a = 9 a (lnp) 
and E Q b = d^ b d a ^(\np) by definition. Following [22[, the 
former describes the density contrast between two neigh- 
boring fluid flow lines, while the latter monitors changes 
(under constant volume) in the shape of an inhomoge- 
neous density distribution (e.g. from spherical symmetry 
to that of an ellipsoidal, etc). Note that, in contrast to 
[lU, here both A a and T, ab are not dimensionless. 

B. Viscous rotating universes 
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Linearized around the Newtonian analogue of the 
Einstcin-de Sitter universe, where matter is in the form 
of pressureless matter (dust), expression (fT4|) reduces to 3 



-2Hu a + vd 2 oj a , 



(16) 



since 0/3 = H in the background. The above implies 
that, within the limits of the linear approximation, vis- 
cosity can no longer generate vorticity, though it has an 
effect on pre-existing rotation. According to (IT"6"1) . in the 
absence of viscous pressure, the vorticity evolution is de- 
termined by the background expansion only. In that case, 
we recover the standard decay-rate of uj a oc a~ 2 oc i~ 4 / 3 , 
since a oc i 2 / 3 in an Einstein-de Sitter universe. The 
introduction of a viscous fluid, however, will generally 
change the standard picture. It is also worth noting that 
expression ([To]) looks like the "heat-transfer equation", 
with an extra term due to the universal expansion and 
the viscosity coefficient (v) playing the role of the thermal 
diffusivity. 

To proceed, we introduce the harmonic splitting uj a = 
^n)Q ( a l \ with d a oj (n) = = Qi n) = d a Q { a n) and 
d 2 Q { a ] = -{n/a) 2 Q ( a l) . Then, Eq. (TjU) recasts into 



J(n) 



-2Huj 



(17) 



where n > is the comoving wavenumber of the vorticity 
mode. Hence, qualitatively speaking, viscosity enhances 
the diluting role of the expansion and leads to less resid- 
ual vorticity. The viscous effect, however, is scale depen- 
dent; see also Eqs. (f2~Tj) and (|22j) . Finally, recalling that 



3 When linearizing, the only zero-order variables are the back- 
ground ones. In an Einstein-de Sitter universe, these are the 
volume scalar (which coincides with the Hubble parameter — 
i.e. H = 0/3), the density and the isotropic pressure of the 
fluid. All these have only temporal dependence. Physical quan- 
tities that vanish in the background, but have nonzero perturbed 
value, like the vorticity, the viscous pressure, the shear and all 
the spatial gradients, are treated as first-order (linear) pertur- 
bations. Note that differentiating perturbations does not affect 
their perturbative order. Finally, products between linear quan- 
tities are of higher-order and are therefore neglected. 



FIG. 1: Evolution of rotational modes, with respect to time, 
on different scales; see Eq. (|19p . Sufficiently small scales are 
strongly affected by viscosity and there vorticity dissipates 
quickly. On larger wavelengths, however, rotational distor- 
tions can survive for long. There, at late times, vorticity is 
only affected by the background expansion. Note that the 
Ay/An = line corresponds to the standard ui oc t" 4/3 evolu- 
tion law, while the A„/A n = 1 curve defines the viscous scale. 



a oc t 2 ! 3 and H — 2/3t in the background - see (fTU|) , we 
arrive at 



W(„) 



3t 



■ UJ 



{n) 



4/3 



UJ 



in) 



(18) 



with the zero suffix indicating a given time during the 
evolution. Clearly, at sufficiently late times the last term 
on the right-hand side of the above becomes negligible 
and the viscosity effects fade away. More specifically, 
assuming that cj(t = t ) = uj , the solution of (fT81 reads 




1- - 



1/3 




(19) 



where \ u = \J vAh /2c is a characteristic length that de- 
termines the domain of viscous influence; see Eq. (|22[) 
below. Also Ah — c/H (with c being the speed of light) 
represents the Hubble horizon and A„ = a/n is the physi- 
cal scale of the perturbation. In solution (fiT)|) we see again 
that at late times the vorticity evolution is determined 
by the background expansion (see Fig. [JJ. In particular, 
when t 3> to, the above expression is approximated by 



uj = ujq exp 



A, 

An/ 



4/3 



(20) 



ensuring that ui oc t -4 / 3 . On the other hand, at suffi- 
ciently early times (i.e. for t <C to), the rotational be- 
havior of the fluid is mainly dictated by the (viscous) 
exponential term of (TT91 . The overall effect, however, 
also depends on the scale of the perturbation. 
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C. The viscous scale 

The scale dependence of the dissipative viscous effects 
is a key feature of solution (TIT))) , because it determines 
their domain of influence and strength. As one can im- 
mediately see in Fig. [T] the smaller the scale of the mode 
the stronger the viscous dissipation. Further understand- 
ing can be achieved through a simple qualitative analysis, 
without even solving Eq. (fT5|). To be precise, recalling 
that Xh = c/H and A„ = a/n, expression (|f 71) recasts as 



-2H 1 



H 



2cXl 



(21) 



Accordingly, the viscous effects are important when 
vXh/2c\\ ^> 1, which holds for modes with wavelength 



X n <C X v 



vX 



H 



2c 



(22) 



The right-hand side of the above defines a characteristic 
length, the "viscous scale" (Ay), beyond which viscos- 
ity is essentially negligible and below which it becomes 
progressively more important. 4 Using simple kinetic- 
theory arguments, one finds that v ~ vXf, where v is 
the random velocity of the particles involved and A/ 
their mean free-path (see [23j j for further discussion and 
details). Hence, X v <C Xh for non-relativistic species, 
which implies that the typical domain where the viscous 
effects can play a significant role is rather small. Fur- 
thermore, following Eq. (f2"2"j) . the viscosity scale grows as 

X v oc X]( 2 oc t 1 / 2 oc a 3 / 4 . On the other hand, A„ oc a 
always. This means that even wavelengths that are ini- 
tially smaller than the viscous length will eventually cross 
outside. 

The effect of fluid viscosity, with respect to the scale 
of the rotational distortion at different times, is plotted 
in Fig. [5] There, one can see that viscosity dictates the 
evolution of rotational perturbations that lie deep inside 
the viscous scale, while much larger lengths are only af- 
fected by the background expansion. To be precise, when 
X u /X n ^> 1, we find that 



wo e x P ' 



1 



1/3 



(23) 



Therefore, modes that have remained well inside the vis- 
cous scale long enough have 



uj = ljq exp 



I [ ± 

A*, 



(24) 



when they cross A„. Afterwards, the background expan- 
sion dominates and these modes evolve according to so- 



lution ([20]) . In other words, the factor exp 



-4(A v /A„)o 



1.0000 




0.0001 



10.0 



FIG. 2: Evolution of vorticity perturbations, with respect to 
the scale of the mode, at different times (see Eq. ()19p ). Ini- 
tially, when t/to ~ 1, all wavelengths have essentially the 
same amount of rotation (uj/ojq ~ 1 almost everywhere). 
As time goes by, the vorticity decays away. Dissipation is 
stronger on small scales (i.e. for A„/A„ S> 1), where the vis- 
cous effects dominate. On large lengths, where the decay is 
mainly due to the expansion, the vorticity survives longer. 



in Eq. (|2TJ|) expresses the total vorticity reduction that is 
caused by the fluid viscosity on small scales, while the 
power-law reflects its subsequent decay due to the uni- 
versal expansion. Clearly, the dissipative effect of viscos- 
ity on modes that remained inside the viscous scale for 
long (i.e. those with X v /X n 3> 1 initially) is very strong, 
leaving no residual vorticity for all practical purposes. 
Finally, we should also note that, in analogy to Silk- 
damping [24[ , viscous dissipation could potentially affect 
scales considerably larger than A M . Investigating whether 
and how this happens, however, goes beyond the scope 
of the present paper. 



IV. MAGNETIZED ROTATION 

A. Magnetically induced vorticity 

Let us now consider the rotational behavior of a self- 
gravitating Newtonian medium in the presence of a mag- 
netic field. The vorticity evolution is still monitored by 
Eq. iQJ, though the shape of the Navier-Stokes formula 
changes. Assuming zero fluid pressure and very high elec- 
trical conductivity (i.e. adopting the ideal MHD approx- 
imation), we have 



pA a 



ibe 



B"curlB c 



1 



d a B 2 + B b d b B a , (25) 



4 On very small scales, we can no longer ignore the nonlinear effects 
and the linear approximation breaks down. 



where the right-hand side provides the Lorentz force 
and we are employing Heaviside-Lorentz units for the 
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Maxwell field. 5 In many applications it helps to split the 
Lorentz force into a pressure and a tension part, with the 
former reflecting the tendency of the magnetic force- lines 
to push each other apart and the latter their "desire" to 
remain straight. 



Taking the "curl" of (f^5|) and then substituting the 
resulting expression into the right-hand side of (g]) , while 
keeping in mind that d a B a = at the ideal-MHD limit, 
we arrive at (see Appendix iBl for more details) 



u a = - \ ®Lo a - e abc d" (lnp) d c B 2 
6 Ap 

+ Sabcd" (In p) B d d {d B c ^ - -1 cui\B b d b (In p) B a 
2p Ap 



— B b d b (lnp) curLB + — cm\B b d {b B a) 

— B b d b cm\B a + a ab iu b . 
2p 



(26) 



According to the above, which provides an irreducible 
decomposition of the vorticity propagation formula, the 
B-field will generally act as a source of rotation, even in 
the case of a homogeneous density distribution. 6 We also 
note that, of the six magnetic source-terms on the right- 
hand side of (f26]). only the first comes from the pressure 
component of the Lorentz force; see Eq. ([25]) . This term 
may be seen as the magnetic analogue of the baroclinic 
stress that emerges in non-barotropic fluids (see and 
also compare to Eq. (fT4"|) here). The rest of the mag- 
netic source terms in (|2"o) are due to the field's tension. 
All these imply that, in the absence of density inhomo- 
geneities, vorticity is generated solely by the magnetic 
tension. 



tioned relation reduces to 7 

dj a = -2Hu) a + i cm\B b d {b B a) - i B b d b cm\B a . (27) 
2p 2p 

Therefore, large-scale magnetic fields can act as sources 
of linear vorticity in cosmological environments [1 31 ] . The 
difference with the nonlinear expression (1261) given before 
is that now the sole source of vorticity is the magnetic 
tension. In other words, it is the negative pressure that 
the £>-ficld exerts along its own direction which generates 
linear rotational perturbations. We should also note that 
the last term in Eq. (|2"T|) will not lead to vorticity genera- 
tion when the "curl" of the field remains invariant along 
its own direction. This happens, for example, when the 
field lines are straight or circular. 

Expressions (|2"6) and (j2"Tf also allow one to look at 
the magnetic effect on pre-existing vorticity. For exam- 
ple, the presence of the field will generally cause pre- 
cession effects by changing the direction of the vorticity 
vector, namely the rotational axis of the fluid. Further 
insight can be obtained by taking the convective deriva- 
tive of ((27)) . Then, keeping in mind that H = —np/2 
and p = —3Hp in the Einstein-de Sitter background - 
see Eqs. (|Hb) and © - gives 



id a = -5HLJ a - KpUJ a . 



(28) 



Note that in deriving the above we have also used the 
1/2-order relations B a = —2HB a , which is the in- 
duction equation, and the associated commutation laws 
(d b B a y = ~3Hd b B a , (curLB Q )' = -3Hcw\B a and 
(<9;,curl£> a )' = —AHd b cur\B a (see also in Appendix [Cl for 
further details). In addition, to eliminate the magnetic 
terms, we have re-employed expression (|27p. Finally, re- 
calling that H = 2/3t and up — 4 /it 2 in an Einstein-de 
Sitter universe (see Eq. (fTU|)). expression ([28]) reads 



10 



■ r 



t 2 UJ a 



(29) 



The latter accepts the solution 



B. Magnetized rotating universes 



Expression (1261) describes the nonlinear rotational be- 
havior of a self-gravitating, Newtonian medium, of zero 
pressure and high electrical conductivity. When lin- 
earized around the (non-magnetized) Newtonian ana- 
logue of the Einstein-de Sitter universe, the aforemen- 



5 One can also obtain l|25j l from the Navier-Stokes equation of an 
imperfect fluid, namely from mill . Indeed, recalling that the 
B-ficld corresponds to a imperfect medium with p = B 2 /6 and 
^ab = n a6 = (B 2 /3)S ab - B a B b , with d a B a = at the ideal- 
MHD limit, it is straightforward to show that expression (IIH 
leads to Eq. {25jl. 

6 As in Eq. II15H before, one can replace the gradient <9 a (lnp) with 
the vector A a , where the latter defines the (non-dimensionlcss) 
density contrast between two neighboring flow lines. 



0J„ OC t 



- 1 oc a- z ' 2 . 



(30) 



since a oc t 2 / 3 to zeroth order. This result shows that 
magnetized vortices decay inversely proportional to time, 
as opposed to their uj a oc i~ 4 / 3 oc a~ 2 evolution in the 
non-magnetic case. The presence of the i?-field has re- 
duced the decay rate of linear vorticity perturbations. 
In other words, a magnetized Einstein-de Sitter uni- 
verse should contain more residual rotation than its non- 
magnetized counterpart. An analogous conclusion has 



7 When linearizing H26H we have treated the magnetic pressure 
(isotropic and anisotropic) as linear perturbation. This ensures 
that both B 2 and Tl ab are first-order variables, while the mag- 
netic vector and its spatial gradients are 1/2-order perturbations. 
Then, the magnetic terms on the right-hand side of l|27| l have 
perturbative order one. 
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also been reached, through a relativistic analysis, in- 
directly - after studying the magnetic effects on rota- 
tional (vortex-like) density perturbations in an almost- 
FRW universe [Hi . 

Unlike viscosity, the linear magnetic effect on vorticity 
is scale independent. In practice, this means that the as- 
sociated domain of influence is essentially determined by 
the coherence length of the -B-field. Also in contrast to 
viscosity, the magnetic presence seems to "favor" both 
the generation and the survival of linear vorticity per- 
turbations. The latter effect, in particular, should be 
directly attributed to the tension properties of the mag- 
netic field lines, namely to the negative pressure the B- 
field exerted along its own direction. 8 



V. DISCUSSION 

Realistic media do not behave like perfect fluids and 
magnetic fields, of various scale and strength, appear to 
be everywhere in the universe. At the same time, cosmic 
rotation still poses a number of fundamental questions 
that remain as yet unanswered. In this respect, it would 
be interesting to see whether viscosity and magnetism 
have something new to say about vorticity in general. 
So far, both agents have been known to act as sources 
of rotation, though their effects on pre-existing vorticity 
are not clear yet. 

The present article adopts the Newtonian approach to 
investigate the rotational behavior of a viscous fluid and 
subsequently that of a magnetized medium. At first, our 
analysis is fully nonlinear and confirms that both viscos- 
ity and magnetism can generate vorticity. We also show 
that the former of these two agents does so through the 
coupling of the shear with certain types of density in- 
homogeneities. Magnetism, on the other hand, does not 
necessarily need an inhomogeneous matter distribution 
to source rotational distortions. The B-field generates 
vorticity by means of its Lorentz force, especially via the 
tension component of the latter. 

Assuming the Newtonian analogue of the Einstein-de 
Sitter universe and then linearizing our equations around 
this background, we looked into the effects of viscosity 
and magnetism on rotational cosmological perturbations. 
We found that the former effects are negative. Viscos- 
ity generally leads to vorticity dissipation, though (for 
non-relativistic species) its role is typically confined to 
relatively small lengths. Nevertheless, further study is 
required to establish the actual range of the dissipativc 



8 When dealing with linear density (scalar) perturbations, only 
the positive magnetic pressure is accounted for. In that case, 
the effect of the field is to reduce the growth rate of the inho- 
mogeneities, or to increase the associated Jeans scale (i.e. the 
region where pressure support prevents these distortions from 
growing) pjj- 



effects. Here, we have taken the first step in this di- 
rection by identifying the characteristic "viscous scale", 
below which viscosity dominates and beyond which the 
universal expansion takes over. 

In contrast, the magnetic presence seems to favor ro- 
tation. More specifically, after linearizing our equa- 
tions around the aforementioned Einstein-de Sitter back- 
ground, we found that the -B-field slows down the stan- 
dard decay rate of rotational cosmological perturbations. 
In particular, vorticity no longer obeys the standard a -2 
decay-law, but drops as a~ 3 / 2 (a is the cosmological scale 
factor). Also, the magnetic effects are scale-independent, 
which means that their domain of influence is decided 
by the coherence length of the £>-field. Analogous re- 
sults were also obtained when studying vortex-like per- 
turbations in the density distribution of a magnetized 
medium. Overall, it appears as though magnetism fa- 
vors both the emergence and the survival of cosmic ro- 
tation. Put another way, magnetized universes should 
contain more residual vorticity than their magnetic-free 
counterparts. 

In closing, we should also note that it would be worth- 
studying the above described effects using a relativistic 
approach. This will enable one to involve highly ener- 
getic particles, look deep into the pre-decoupling uni- 
verse and also allow for spacetimes more complex than 
the Einstein-de Sitter universe. 
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Appendix A: The viscous source term 

Viscosity generates vorticity through the third term in 
the right-hand side of Eq. (fH| . Further insight can be ob- 
tained by decomposing this source term to its irreducible 
components. To begin with, we have 

e abc d b [<r cd d d (\np)] = s abc d b a cd d d (lnp) 
+e abc a cd d b d d (\np) 
= e abc d b a cd A d 

+e abc o c d Y, db , (Al) 

where A a = 9 a (lnp) and Yi ab = dn,d a \(hip) by defini- 
tion. The antisymmetric part of d b d a (In p) is identically 
zero, while its trace has no contribution to decomposi- 
tion (lAip . due to the total antisymmetry of the Levi- 
Civita tensor and the symmetry properties of the shear. 
Moreover, recalling that cvx\a ab = e cd ( a d c a d b ) and the 
trace-free nature of the shear tensor, we may write 

curkr a6 = i e cda d c a d b + i e cdb d c a d a . (A2) 
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Also, given that curla a b = —d{ b u) a ) (see constraint (JSJa)), 
the above recasts into the auxiliary relation 

A b e cda d c a\ = -2A b d {b iu a) - A b e bcd d c a d a . (A3) 

Finally, substituting this result into the right-hand side 
of (IA1[) and employing some straightforward tensor alge- 
bra, one arrives at decomposition (|15|) . The latter shows 
that only part of the viscous source term in Eq. (|14j) can 
actually generate vorticity. 

Appendix B: Nonlinear magnetized vorticity 

The rotational evolution of a medium is governed by 
Eq. ((4]), where the key term is the "curl" of the accel- 
eration vector. By definition, the latter is curb4 a = 
e abc d b A c . Then, starting from the magnetized Navier- 
Stokes equation (see expression ([25]) ). we first arrive at 

curL4 Q = ±- s abc d b (In p) d c B 2 - - e abc 8 b (In p) B d d d B c 
2p p 

+- e ab c 8 b B d d d B c + - B b d b cur\B a . (Bl) 
P P 

Recalling that d a B a — at the ideal-MHD limit and 
using the total antisymmetry of the Levi-Civita tensor, 
the third term on the right-hand side of the above gives 

e ab c d b B d d d B c = 2e ab c d( b B^d [d B c] 

= -cnx\B b d {b B a) , (B2) 

since that 2d[ a B b ] = £ a ;, c curLB c . In addition, given that 
d b B a = d( b B a ) +d[ b B a ] , the second term in the right-hand 
side of IjBip decomposes to 

e abc d b (\np)B d d d B c = e abc d b (hip) B d d< d B c ^ 
+ \B b d b (lnp) curLB Q 

-i C urLB 6 c> b (lnp) B a .(B3) 

Substituting the last two auxiliary relations back into 
Eq. (jBll) and the resulting expression into (@| we obtain 



Eq. (|26l) . The latter monitors the nonlinear evolution of 
rotational perturbations in a magnetized, highly conduc- 
tive medium. 

Appendix C: Linear commutation laws 

In Euclidean spaces, simple partial derivatives (tem- 
poral or spatial) commute. This is not the case, how- 
ever, when convective derivatives are involved. Then, 
one needs to employ the associated commutation laws. 
In § IIVB( we were only interested in the linear version 
of these formulae. Hence, starting from the convective 
derivative of d b B a , using decomposition (|T|) and then lin- 
earizing around an Einstein-de Sitter universe, one finds 

(d b B a y = d t d b B a + v c 8 c d b B a 

= d b d t B a + d b {v c d c B a ) - d b v c d c B a 

= d b B a - Q QS c b + a% + e c bd uj d ^j d c B a 
= -3Hd b B a . (CI) 

Also note that B a = —2HB a to first order, while ~ 
3H and a ab = = ui a in the unperturbed background. 
Similarly, we have 

(curLB a )' = e abc d t d b B c + e abc v d d d d b B c 

= e abc d b d t B c + e abc d b (v d d d B c ) 

-e abc d b v d d d B c 
= e abc d b B c - HcurlB a 
= -3Hcm\B a , (C2) 

to linear order. Finally, proceeding in an analogous way 
and then combining the last two results, one arrives at 

(d b cm\B a y = -4Hd b mr\B a . (C3) 

On using these commutation laws, between convec- 
tive and spatial partial derivatives, one can begin from 
Eq. ([27]) and obtain expression (|28[) . The latter monitors 
the linear evolution of magnetized rotational perturba- 
tions on an Newtonian Einstein-de Sitter universe. 
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